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EXECUTIVE  SUMMARY 


The  main  objective  of  this  basic  research  program  was  to  develop  a  generic  physics-based 
computational  methodology  for  accurately  predict  the  fatigue  properties  of  structural  alloys  and  their 
variability  with  microstructure.  The  objective  of  the  program  was  achieved  in  three  tasks  with  these 
aims:  (1)  to  develop  physically  based  fatigue  crack  initiation  and  growth  models,  (2)  to  develop  a 
probabilistic  approach  for  linking  physically  based  models  into  a  continuum  framework,  and  (3)  to 
demonstrate  the  utilities  of  the  methodology  in  a  probabilistic  design  setting 

An  integrated  interdisciplinary  modeling/experimental  approach  involving  a  combination  of 
materials  sciences,  micromechanics,  and  probabilistic  techniques  was  used  to  develop  a  suite  of 
physically-based  models  for  describing  dislocation  structure  evolution,  fatigue  crack  initiation, 
fatigue  crack  growth  (FCG)  and  threshold  behaviors,  as  well  as  the  linkage  of  these  processes  over 
the  atomistic,  dislocation  cell,  grain,  and  continuum  size  scales.  This  set  of  fatigue  models  was 
incorporated  into  a  probabilistic  framework  and  integrated  into  a  computer  code  named 
MicroFaVa  (Micromechanical  Fatigue  Variability  Code)  for  predicting  fatigue  properties  and  life 
variability  due  to  material  variation.  The  probabilistic  material  code  was  coupled  to  and  utilized 
with  DARWIN®,  a  probabilistic  design  and  life-prediction  code  developed  at  Southwest  Research 
Institute  (SwRI®),  to  treat  fatigue  life  variability  in  a  Ti  rotor  due  to  variation  in  the  microstructure. 

The  key  accomplishments  of  this  program  are  (1)  the  development  of  physically-based 
fatigue  crack  initiation  and  growth  models;  (2)  the  development  of  a  probabilistic  code  named 
MicroFaVa  for  treating  fatigue  life  variability  resulting  from  material  variation;  (3)  the  incorporation 
of  MiCroFaVa  with  a  probabilistic  design  code  DARWIN®  and  their  application  in  the  life  prediction 
for  a  structural  component  made  of  Ti-6A1-4V;  (4)  an  extensive  and  systematic  assessment  of  the 
predictive  capability  of  the  physics-based  fatigue  models;  and  (5)  verification  and  validation  of  the 
integrated  probabilistic  design  and  lifing  methodology  by  applying  the  computational  approach  to 
two  benchmark  problems  with  known  solutions.  In  addition,  the  utility  of  the  computed  material 
input  was  successfully  demonstrated  by  one  to  one  comparisons  against  the  conventional  method 
based  on  empirical  data  input.  The  probabilistic  micromechanical  code  was  utilized  to  elucidate  the 
influence  of  microstructural  variation  on  the  propagation  and  arrest  of  grain-sized  small  cracks,  the 
probability  density  function  of  small-crack  thresholds,  as  well  as  the  occurrence  of  dual  fatigue  limits 
for  crack  initiation  and  growth.  These  computational  results  were  used  to  develop  fundamental 
understandings  on  pertinent  microstructure/fatigue  resistance  relationships  such  as:  (i)  the 
dependence  of  the  intrinsic  threshold  on  the  planarity  of  slip  and  dislocation  cell  formation,  (iii)  the 
occurrence  of  dual  fatigue  thresholds  and  dual  fatigue  limits  resulting  from  differences  in  fatigue 
mechanisms,  (iii)  fatigue  life  variation  resulting  from  difference  in  the  slip  morphology  and  enhanced 
fatigue  crack  initiation  life  by  dislocation  cell  formation,  and  (iv)  dependence  of  FCG  rate  on 
dislocation  density  at  the  crack  tip.  Most  of  these  results  are  published  in  the  open  literature  or  have 
been  submitted  for  publication. 

The  technology  developed  in  this  program  has  been  incorporated  into  DARWIN®  and 
transferred  to  the  USAF  DUST  program  (Contract  No.  F336 15-03-2-5203)  at  SwRI.  The 
probabilistic  micromechanics-based  life  methodology  is  in  the  process  of  being  transferred  to 
a  DARPA-supported  program  entitled  “A  Revolutionary  Smart  Materials  Approach  to  Structural 
Health  Monitoring  and  Progress”  (Contract  No.  HR001 1-04-0004)  at  SwRI  and  an  AFOSR  program 
on  Vehicle  Health  Monitoring  (Contract  No.  FA955004 10254)  at  University  of  Texas  at  San 
Antonio. 
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1.0  INTRODUCTION 


There  is  a  need  to  develop  affordable  new  materials  by  reducing  the  long  lead  time  and 
the  number  of  material  and  processing  iterations  required  to  achieve  the  desired  material 
properties.  Traditional  approaches  based  on  trails  and  errors  used  in  material  development  need 
to  be  replaced  with  computer-assisted,  science-based  approaches  so  the  desired  mechanical 
properties  can  be  achieved  in  new  materials  systematically,  efficiently,  and  affordably.  Thus, 
there  is  a  need  for  a  generic  physics-based  computational  methodology  that  accurately  predicts 
the  mechanical  properties  of  new  materials  and  their  microstructure  so  that  computer-assisted 
optimization  procedure  can  be  applied  to  obtain  the  desired  microstructures  and  mechanical 
properties. 

Fatigue  crack  growth  properties  of  structural  alloys  generally  exhibit  variability  because 
of  variations  in  microstructure.  This  variability  is  traditionally  characterized  by  extensive  testing 
to  develop  a  large  database  so  that  the  mean,  variance,  and  the  lower  bound  of  the  fatigue 
properties  can  be  ascertained.  An  obvious  drawback  of  this  approach  is  that  an  experimental 
database  is  time-consuming  and  costly  to  generate.  Another  limitation  is  that  the  empirical 
approach  seldom  provides  information  regarding  the  sensitivities  of  various  material  parameters 
to  the  probability  of  fracture  associated  with  a  particular  component  or  the  means  for  reducing 
fatigue  variability.  One  of  the  contributing  factors  to  this  problem  is  that  most  existing  fatigue 
life  or  crack  growth  models  are  not  linked  with  the  fundamental  material  parameters;  instead 
they  are  based  on  empirical  constants  obtained  from  curve  fits  of  data  from  physical  tests. 
Although  these  models  characterize  the  main  descriptors  (mean,  variance)  of  fatigue  life  as  a 
function  of  applied  stress  range  (S-Nf  curves)  or  fatigue  crack  growth  rate,  da/dN,  as  a  function 
of  stress  intensity  factor  range,  AK,  they  provide  little  or  no  physical  insight  into  the  material 
parameters  that  influence  failure  associated  with  fatigue  crack  nucleation  and  growth. 

In  general,  fatigue  life  can  generally  be  considered  to  consist  of  two  parts,  which  are 
crack  initiation  and  crack  propagation  lives.  The  crack  propagation  life,  Nf,  can  be  predicted  by 
integrating  the  da/dN  equation,  leading  to 

*  ■  r; 

where  a,  and  a/are  the  initial  and  final  crack  length,  respectively,  and  a/  is  defined  as  the  critical 
crack  length  at  fracture.  The  fatigue  crack  growth  (FCG)  life,  Nf,  is  most  sensitive  to  the  initial 
crack  size,  a,  [1],  but  also  depends  on  microstructure,  which  can  alter  intrinsic  material  properties 
as  well  as  extrinsic  properties  originating  from  crack  shielding  mechanisms. 

As  a  result,  a  physics-based  computational  methodology  required  to  address  the  effects 
of  microstructure  variation  on  fatigue  life  variability  must  include:  (1)  a  microstructure-based 
crack  initiation  model  for  predicting  initiation  life  and  initial  crack  size  distribution,  (2)  a 
microstructure-based  fatigue  crack  growth  model,  and  (3)  a  probabilistic  framework  for  treating 
random  material  variables  in  the  fatigue  models  and  for  probabilistic  design  and  life  assessment 
of  structural  components.  This  report  describes  the  result  of  a  basic  research  program  to  develop 
a  generic  physics-based  computational  methodology  for  predicting  fatigue  life  response  of 
structural  alloys  and  their  variability  with  microstructure. 
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2.0  PROGRAM  OBJECTIVES  AND  APPROACH 
2.1  Research  Objectives 

The  main  objective  of  this  basic  research  program  was  to  develop  a  generic  physics-based 
computational  methodology  for  accurately  predict  the  fatigue  properties  of  structural  alloys  and  their 
variability  with  microstructure.  The  objective  of  the  program  was  achieved  in  three  tasks  with  these 
aims:  (1)  to  develop  physically  based  fatigue  crack  initiation  and  growth  models,  (2)  to  develop 
a  probabilistic  approach  for  linking  physically  based  models  into  a  continuum  framework  and  a 
generic  physics-based  computational  code  for  predicting  fatigue  life  properties  and  variability  in 
structural  alloys,  and  (3)  to  demonstrate  the  utilities  of  the  methodology  in  a  probabilistic  design 
setting. 


2.2  Approach 

An  integrated  interdisciplinary  modeling/experimental  approach  involving  a  combination 
of  materials  sciences,  micromechanics,  and  probabilistic  techniques  was  used  to  develop  a  suite 
of  physically-based  models  for  describing  dislocation  structure  evolution,  fatigue  crack  initiation, 
fatigue  crack  growth  (FCG)  and  threshold  behaviors,  as  well  as  the  linkage  of  these  processes 
over  the  atomistic,  dislocation  cell,  grain,  and  continuum  size  scales. 


3.0  SUMMARY  OF  KEY  RESULTS 

The  key  results  of  this  program  are  (1)  the  development  of  physically-based  fatigue  crack 
initiation  and  growth  models,  (2)  the  development  of  a  probabilistic  code  named  MicroFaVa  for 
treating  fatigue  life  variability  resulting  from  material  variation,  and  (3)  the  incorporation  of 
MiCroFaVa  with  a  probabilistic  design  code  DARWIN  and  their  application  in  the  life 
prediction  for  a  structural  component  made  of  Ti-6A1-4V.  In  this  report,  a  brief  summary  of  the 
probabilistic  fatigue  crack  initiation  and  growth  models  in  MicroFaVa  is  first  presented.  Second, 
the  application  of  the  probabilistic  model  to  predicting  the  FCG  variability  in  Ti-6A1-4V  is 
highlighted.  Finally,  application  of  the  model  to  lifing  a  Ti  rotor  design  is  demonstrated. 


3.1  Probabilistic  Fatigue  Crack  Growth  Rate  Modeling 

The  physics-based  crack  growth  model  proposed  by  Chan  and  Enright  for  treating  the 
three  stages  of  fatigue  crack  growth  has  the  form  given  by  [5,  6,  7]: 
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where  the  fatigue  crack  growth  rate,  da/dN,  is  expressed  in  terms  of  an  effective  stress  intensity 
range,  AKeff,  which  is  the  difference  between  the  applied  stress  intensity  range,  AK,  and  the  stress 
intensity  range  at  crack  closure,  AKci\  £  a  dimensionless  normalizing  parameter  defined  in  terms 
of  the  Young’s  modulus  (£),  dislocation  cell  size  (s),  cyclic  yield  stress  (<Jy),  fatigue  ductility 
coefficient  (£ /),  and  the  dislocation  barrier  spacing  (d).  The  dislocation  barrier  spacing  is  taken 
to  be  a  function  of  the  grain  size  D  as  described  byEq.  (4),  where  d„,  D,„  and  y  are  material 
constants.  The  stress  intensity  range  at  the  Stage  I  to  Stage  II  transition,  AKt,  is  a  function  of  the 
fatigue  limit  (Acre),  cyclic  yield  stress  range  and  the  large-crack  FCG  threshold  ( AK,h)',  n /  and  «2 
are  the  Stage  I  and  Stage  II  exponents,  respectively.  R  is  the  stress  ratio  and  kc  is  the  fracture 
toughness. 


The  effective  stress  intensity  range,  AKejf,  is  expressed  as  [9] 

A  Krff  =  UfUdUkU0AK  (6) 

where  U,  is  ratio  of  AKej?  to  AK  for  individual  crack  closure  mechanisms  including  plasticity-  ( P ), 
deflection-  (£>),  roughness-  (/?),  and  oxide-induced  (0)  closure.  Analytical  expressions  for 
computing  individual  (/,  terms,  which  are  summarized  in  a  recent  paper  [9],  include  Newman’s 
formulation  for  treating  plasticity-induced  crack  closure  [10],  Suresh’s  formulation  for  crack 
deflection  and  roughness-induced  crack  closure  [1 1],  as  well  as  an  oxide-wedging  model. 

3.2  Probabilistic  Crack  Initiation  Modeling 


The  crack  initiation  model  of  Tanaka  and  Mura  [2]  was  developed  based  on  a  dislocation 
dipole  mechanism  operating  in  a  surface  grain.  During  fatigue  loading,  irreversible  slip  occurs 
on  parallel  slip  planes  in  a  favorably  oriented  surface  grain,  producing  dislocation  dipoles  at  the 
ends  of  a  double  pileup  whose  coalescence  ultimately  leads  to  crack  nucleation.  Chan  [3] 
recently  extended  Tanaka  and  Mura’s  model  [2]  to  explicitly  incorporate  the  crack  size,  as  well 
as  other  pertinent  material  parameters  in  the  response  equation  by  considering  the  energetics  of 
the  fatigue  crack  initiation  process.  Specifically,  the  length  of  the  incipient  crack  was  obtained 
by  equating  the  elastic  strain  energy  released  by  dislocation  coalescence  and  crack  opening  to  the 
fracture  energy,  consisting  of  elastic  and  plastic  components,  required  to  form  the  crack  surfaces. 
This  formulation  leads  one  to  [3,  6] 


{Aa-lMk)N° 


8 M\/d  TVftYaV 

/br(l-v)J  IdAoJ 


(7) 


for  crack  initiation  at  slip  bands,  where  Ac  is  the  stress  range,  M  is  the  Taylor  factor,  k  is 
the  friction  stress,  N-,  is  the  cycles-to-initiation,  v  is  Poisson’s  ratio,  p  is  shear  modulus, 
and  A  (=  0.005)  is  a  universal  constant.  Eq.  (7)  relates  the  crack  initiation  life,  to  the  grain 
size,  D,  the  slipband  width,  h,  and  the  crack  half-length  or  depth,  a.  The  exponent  to  Nj  has  been 
generalized  to  or  where  0  <  a<  1.  The  value  of  oris  not  a  constant,  but  depends  on  the  degree  of 
slip  irreversibility  and  the  stacking  fault  energy.  Expressions  for  crack  initiation  at  inclusions 
and  notches  have  also  been  developed  and  reported  earlier  [3]. 
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3.3  Probabilistic  Framework 


For  probabilistic  modeling,  the  material-specific  variables  in  Eqs.  2  and  7  are  treated 
separately  from  the  remaining  variables.  Previous  studies  (Chan  and  Torng[l],  Enright  and 
Chan  [7])  suggest  that  some  of  the  material-specific  variables  can  be  approximated  as 
deterministic  variables  (e.g.,  E,  s ,  Kc).  The  remaining  material  variables  (dislocation  barrier 
spacing  d,  cyclic  yield  stress  a'y,  fatigue  ductility  coefficient  e'j)  are  modeled  as  random 
variables.  In  many  instances,  these  material  variables  are  all  functions  of  the  grain  size  D.  If  the 
non-material-specific  random  variables  (a„  AK,  R )  are  temporarily  treated  as  deterministic 
variables,  Eq.  (1)  can  be  used  to  predict  Nf  scatter  attributed  to  material  variability.  A  general 
probabilistic  framework  has  been  developed  to  address  the  influences  of  both  the  material-  and 
non-material-specific  variables  on  the  probability  of  fatigue  fracture  (Enright  and  Wu  [12], 
Leverant  [13],  Wu  et  al.  [14],  Enright  et  al.  [15]).  The  Nf  scatter  results  shown  in  this  report  can 
be  directly  applied  to  this  framework,  allowing  the  designer  to  quantify  the  influences  of  the 
fundamental  micromechanics  variables  on  overall  component  risk.  Figure  1  shows  the 
schematics  of  the  probabilistic  framework  utilized  in  MicroFaVa  for  modeling  random  variables. 


Figure  1.  Schematics  of  the  probabilistic  framework  in  MicroFaVa  for  modeling 
material-specific  random  variables. 
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As  in  probabilistic  crack  growth  modeling,  material  variables  in  Eq.  (7)  are  treated 
separately  from  the  remaining  (e.g.,  loading)  variables  in  probabilistic  crack  initiation  modeling. 
The  material  specific  variables  include  grain  size  (D),  slipband  width  (h),  fatigue  limit  (Mk), 
elastic  properties,  and  Taylor  factor,  which  is  influenced  by  texture.  The  remaining  variables  are 
the  stress  range,  cycles-to-crack  initiation  (Nj),  and  crack  size  (a,)  at  initiation.  Eq.  (7)  can  be 
utilized  in  two  ways:  (1)  for  computing  the  Nj  distribution  for  a  given  value  of  crack  size  at 
initiation,  a„  and  (2)  for  computing  the  distribution  of  a,  for  a  given  number  of  fatigue  cycles. 
For  both  cases,  the  general  probabilistic  framework  described  by  Enright  et  al.  [15]  is  employed. 

3.4  Model  Application  to  Ti-6AI-4V 

The  probabilistic  microstructure-based  fatigue  crack  initiation  and  growth  modeling 
methodology  is  illustrated  for  Ti-6A1-4V.  Figure  2  shows  the  typical  microstructure  of 
the  Ti-6A1-4V  alloy,  which  is  comprised  of  primary  or  grains  and  Widmanstatten  a+  (5  colonies. 
The  volume  fraction  of  gc  grain  was  about  60%  and  the  average  grain  size  was  about  12  |im;  both 
were  determined  using  metallographic  techniques.  In  grain  size  measurement,  both  the  a  grain 
size  and  the  a- v  (5  colony  size  were  treated  as  equivalent  and  fitted  to  one  single  probability 
density  function.  A  log-normal  distribution  was  used  to  describe  the  grain  size  distribution  (D) 
ofTi-6Al-4V.  As  shown  in  Figure  3,  this  distribution  provides  conservative  values  of  D 
compared  to  experimental  values. 


Figure  2.  Microstructure  of  Ti-6AI-4V  shows  60%  primary  a  grain  (light  phase) 
and  40%  of  a+  ft  Widmanstatten  colonies. 
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Material  constants  in  the  FCG  model  were  determined  from  FCG  data  for  R  =  0.8  to 
ensure  that  the  FCG  response  was  free  from  crack  closure  effects.  Values  for  the 
material-related  variables  are  indicated  in  Table  1.  Once  intrinsic  material  properties  were 
obtained,  the  FCG  model  was  utilized  to  predict  deterministic  da/dN  curves  for  various  R  ratios 
by  applying  the  crack  closure  model  and  using  average  values  of  grain  size,  yield  stress,  and 
fatigue  ductility  coefficient.  The  crack  closure  model  [9]  includes  plasticity  [10],  deflection-, 
and  roughness-induced  [11]  closure  mechanisms,  but  not  oxide-induced  crack  closure.  For 
crack  growth  in  Ti-6A1-4V,  grain  size  was  the  dominant  microstructural  variable,  and  was 
modeled  as  a  lognormal  random  variable  (mean  =  1 1.7  pm,  standard  deviation  =  3.2  pm)  based 
on  grain  size  measurements  shown  in  Figure  3.  Cyclic  yield  stress  a'y  and  fatigue  ductility 
coefficient  e'f  were  modeled  as  deterministic  variables  to  emphasize  the  influence  of  grain  size 
variability  on  da/dN  and  Nf  scatter  results. 
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Figure  3.  Grain  size  of  Ti-6AI-4V  compared  against  the  log-normal  distribution 
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Table  I.  Values  of  Deterministic  Material  Variables  forTi-6AI-4V  at  24°C  used  in 
Microstructure-based  Fatigue  Crack  Growth  Model  (Eq.  (2)) 


Variable 

Description 

Units 

Value 

E 

Young’s  modulus 

MPa 

1.16x10  s 

ni 

Stage  I  exponent 

— 

23.0 

n2 

Stage  II  exponent 

— 

3.87 

s 

dislocation  cell  size 

fJm 

0.25 

do 

reference  dislocation  barrier  spacing 

jJm 

1.0 

Do 

reference  grain  size 

/ im 

1.0 

R 

stress  ratio 

— 

0.1 

Kc 

fracture  toughness 

MPa  Vm 

66.7 

AKth 

fatigue  threshold 

MPa^Jm 

2.0 

D 

grain  size 

jUm 

11.7 

(J'y 

cyclic  yield  stress 

MPa 

909.33 

£{ 

fatigue  ductility  coefficient 

— 

0.0389 

y 

dislocation  barrier  spacing  exponent 

— 

1.0 

3.4.1  S-Nf  Variability 

The  microstructure-based  fatigue  crack  initiation  and  growth  models  were  used  to  predict 
the  fatigue  cycles  for  initiating  a  thumb-nail  crack  of  length  2a  and  depth  c  and  the  number  of 
cycles  for  this  crack  to  reach  a  critical  depth.  Values  of  material  constant  in  the  crack  initiation 
model  are  presented  in  Table  II.  The  computed  crack  initiation  life,  Nt,  and  crack  growth 
life,  Ng,  for  various  alternating  stresses,  oa,  are  compared  against  experimental  data  in  Figure  4. 
The  corresponding  crack  length  versus  fatigue  cycle  curves  of  crack  initiation  (Nj)  and  crack 
growth  (Ng )  for  a  stress  range  of  552  MPa  at  a  stress  ratio  R  of  0. 1  are  presented  in  Figure  5.  For 
both  cases,  the  fatigue  crack  growth  lives,  Ng,  were  computed  without  a  large-crack  threshold. 
Variations  of  the  Ng  value  were  obtained  when  a  fatigue  crack  growth  threshold  was  employed. 
The  life  variation  depends  on  the  crack  size  at  which  a  crack  is  considered  to  have  “initiated”. 
Different  curves  of  crack  length  versus  fatigue  cycle  were  predicted  for  crack  initiation  and 
growth.  For  most  cases,  the  transition  from  crack  initiation  to  crack  growth  appeared  to  occur  at 
a  crack  size  of  one  to  two  grain  diameters,  as  shown  in  Figures  4  and  5. 
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Table  II.  Values  of  Deterministic  Material  Variables  for  Ti-AI-4v  at  24°C  used  in 
Microstructure-Based  Fatigue  Crack  Initiation  Model  (Eq.  (7)) 


Variable 

Description 

Units 

Value 

M 

Taylor  Factor 

— 

2 

A 

Shear  Modulus 

MPa 

4.4x1 04 

V 

Poisson’s  Ratio 

— 

0.333 

A 

Universal  Constant 

— 

0.005 

a 

Fatigue  Initiation  Life  Exponent 

— 

0.5 

oe 

Fatigue  Limit  (Mk) 

MPa 

200  -  272.5 

H 

Slipband  Width 

m 

5x1 0'8 

The  microstructure-based  fatigue  crack  initiation  and  growth  models  were  utilized  to 
predict  the  variability  of  fatigue  crack  initiation  and  growth  live  due  to  grain  size  variations.  A 
comparison  of  the  predicted  and  observed  probability  density  functions  (PDF)  for  crack  initiation 
life  is  shown  in  Figure  6,  which  shows  that  the  range  of  initiation  life  values  associated  with  the 
experimental  data  is  within  the  values  predicted  by  MicroFaVa. 


Figure  4.  Predicted  crack  initiation  life,  A/,,  crack  growth  life,  Ng,  and  total  life,  Nf, 
compared  against  observed  failure  life  for  Ti-6AI-4V.  Experimental  data 
are  from  [17].  From  Chan  et  al.  [16]. 
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Fatigue  Cycles 


Figure  5.  A  comparison  of  predicted  and  measured  crack  lengths  versus  fatigue 
cycles  for  the  initiation  (Nl)  of  one  grain-sized  crack  and  its  growth  ( Ng )  to 
failure.  Experimental  data  are  from  [17].  From  Chan  et  al.  [16]. 


3.4.2  Fatigue  Threshold  Variability 

The  FCG  threshold  model  in  MicroFaVa  contains  an  intrinsic  term  which  originates  from 
the  crack-tip  cyclic  slip  process,  and  an  extrinsic  term  which  arises  from  a  combination  of 
plasticity-induced  closure  with  crack  deflection,  asperity-induced,  and  oxide-induced  crack 
closure  mechanisms.  In  particular,  the  AK,h  expression  is  given  by  [9] 

AK,„  =U-'  AKlhM  =  [U„  U„  UR  U0Y  A K,hJn  (8) 

with  the  intrinsic  threshold,  AKth.m  given  by  [9] 


Ebv 2 

where  b  is  the  magnitude  of  the  Burgers  vector.  The  proposed  model  was  applied  to  compute 
the  AK,h  for  Ti-6A1-4V  for  R  values  ranging  from  -1  to  1  and  the  results  are  presented  in 
Figure  7.  For  these  calculations,  AK,h,m  ranges  from  1 .2  MPaVm  to  2.0  MPaVm.  Most  of  the 
crack  closure  was  contributed  by  crack-wake  plasticity.  Figure  7  shows  a  comparison  of  model 
calculations  against  experimental  data  of  Ti-6A1-4V  from  the  HCF  program  [  1 8]  at  Air  Force 
Research  Laboratory  (AFRL),  Marci  et  al.  [19,  20],  and  Boyce  and  Ritchie  [21]. 
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Cycles  to  Initiation,  x  106 

Figure  6.  Measured  PDF  for  cycles-to-initiation  compared  against  MicroFaVa 
model  prediction  based  on  measured  grain-size  distribution.  From 
Chan  et  al.  [16]. 
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Comparison  of  model  predictions  against  measured  AK,h  for  Ti-6AI-4V  as 
a  function  of  Kmax.  From  Chan  [9]. 


TheTi-6AI-4V  alloy  studied  by  Boyce  and  Ritchie  [21]  were  identical  to  those  used  in 
the  HCF  program;  hence,  there  was  good  agreement  between  these  two  sets  of  AK,i,  data.  For 
these  two  sets  of  data,  AK,ihin  =  2  MPaVm,  it  was  correctly  predicted  by  the  intrinsic  threshold 
model,  Eq.  (9).  The  increase  in  AK,h  at  R  <  0.8  was  caused  entirely  by  crack  closure 
mechanisms.  The  lower  AK^n  value  (1.2  MPaVm)  for  the  Ti-6A1-4V  material  studied  by  Marci 
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et  al.  [19,  20]  was  predicted  from  Eq.  (9)  using  a  Taylor  factor  of  M  =  2.  This  value  of  M 
appeared  to  give  the  lower  bound  for  AKlhj„ .  Again,  the  increase  in  AK,h  at  R  <  0.5  was  entirely 
due  to  crack  closure  mechanisms.  Thus,  the  variation  of  AK,h,in  from  1.2  MPaVm  to  2  MPaVm 
could  arise  from  microstructural  effects  such  as  texture  (different  values  of  M),  yield  stress,  and 
cell  formation.  On  the  other  hand,  the  variation  of  AK,h  at  low  R  ratio  ( R  <  0.5)  is  likely  caused 
by  crack  closure. 

3.4.3  Crack  Growth  Rate  Variability 

Crack  growth  rate  (mean  and  l%/99%  confidence  limit)  values  for  R  =  0.1  predicted 
using  Eq.  (2)  are  compared  against  experimental  data  [21-24]  in  Figure  8.  A  total  of  67  sets 
of  da/dN  data  of  Ti-6A1-4V  of  the  same  or  equivalent  material  and  microstructure  are  presented 
in  Figure  8.  Experimental  da/dN  data  at  selected  A K  values  were  also  analyzed  to  obtain 
probability  densities.  These  results  are  compared  against  probability  densities  and  coefficient  of 
variation  (COV)  of  predicted  crack  growth  rate  at  selected  AK  values  in  Figure  9.  It  can  be 
observed  from  Figure  9(a)  that  the  PDF  values  based  on  the  experimental  data  are  bound  by  the 
microstructure-based  PDF  values  for  most  of  the  AK  values  considered.  Figure  9(b)  shows  that 
the  predicted  COV  is  about  0.54  -  0.56  over  the  AK  range  of  3  MPaVm  to  90  MPaVm.  The  two 
largest  values  of  the  experimentally  determined  COV  are  0.45  at  AK  =  15  MPaVm  and  0.48 
at  AK  =  45  MPaVm.  The  larger  discrepancy  at  AK  values  in  the  range  of  18-30  MPaVm  is  due 
to  insufficient  experimental  data  in  that  crack  growth  regime. 
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Figure  8.  Measured  da/dN  data  compared  against  model  predictions  based  on  a 
log-normal  grain  size  distribution  [16]. 
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Figure  9. 


AK,  MPa  V m 

(a) 


(b) 

Comparison  of  computed  crack  growth  rate  predictions,  probability  density 
functions,  and  coefficient  of  variation  (COV)  with  experimental  data 
for/?  =0.1  at24°C:  (a)  da/dN  and  probability  density  functions,  and  (b) 
coefficient  of  variation  of  da/dN. 
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3.4.4  Crack  Growth  Life  Variability 

Computational  and  experimental  da/dN  values  were  applied  to  the  life  prediction  of 
a  10  cm  x  10  cm  Ti-6A1-4V  plate  with  material  properties  indicated  in  Table  1.  A  uniform 
stress  range  of  600  MPa  was  applied  to  the  plate  at  a  stress  ratio,  /?,  of  0. 1  and  an  initial 
crack  (0.035  cm  x  0.035  cm)  was  placed  at  one  of  the  comers.  Crack  growth  life  was  computed 
using  probabilistic  fracture  mechanics  software  (DARWIN®  2004  [8])  using:  (1)  tabular  da/dN 
values  based  on  Eq.  (2),  and  (2)  da/dN  slope  and  intercept  constants  based  on  linear  and 
bilinear  curve  fits  of  experimental  data  from  the  literature  (Raizenne  1987  [22],  Mom 
and  Raizenne  1988  [23]).  The  results  are  shown  in  Figure  10.  It  can  be  observed  that  all  of  the 
predicted  crack  growth  life  values  based  on  the  experimental  da/dN  data  fall  within  the  5% 
and  95%  confidence  limits  based  on  life  values  associated  with  the  computational 
model  (Eq.  (2)).  The  corresponding  computed  and  observed  probability  density  of  crack 
propagation  life  are  compared  in  Figure  11,  which  shows  that  the  predicted  mean  is  in  good 
agreement  with  the  experimental  data.  On  the  other  hand,  the  predicted  variability  is  larger  than 
experimental  observations  and  is,  thus,  more  conservative.  As  shown  in  Figure  1 1 ,  the  observed 
COV  is  about  0.23  while  the  predicted  COV  is  0.57.  The  observed  COV  is  lower  partly  because 
of  the  sampling  size  and  partly  because  the  COVs  associated  with  curve-fitting  of  individual 
da/dN  datasets  have  not  been  included  in  the  computation  of  the  experimental  COV.  The 
observed  COV  is  expected  to  approach  the  predicted  COV  as  the  sampling  size  increases  and  the 
COVs  of  individual  da/dN  datasets  are  taken  into  account.  Additional  FCG  life  computations 
for  A  a  =  552  MPa  indicates  that  the  predicted  COV  is  independent  of  the  applied  stress  and 
depends  only  on  the  grain  size  distribution  since  stress  variations  are  not  considered  in  these 
calculations. 


Figure  10.  Comparison  of  predicted  crack  growth  life  values  based  on 
microstructure-based  computational  model  and  experimental  data  da/dN 
results  [22,  23]  for  R  =  0.1  at  24°C  [16]. 
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Figure  1 1 .  Comparison  of  computed  and  observed  probability  density  of  FCG  lives 
based  on  microstructure-based  computational  model  and  experimental 
data  da/dN  results  [22,  23]  for  R  =  0.1  at  24°C. 


3.5  Application  to  Probabilistic  Life  Prediction  of  Titanium  Components 

Two  example  problems  were  solved  to  validate  the  MicroFaVa  prediction  and  its 
integration  with  DARWIN®.  For  the  first  validation  case,  the  calibration  problem  of  a  fictitious 
titanium  rotor  disk  subjected  to  a  single  inertial  load  in  FAA  Advisory  Circular  AC33.14-1  [25] 
and  an  empirically  derived  hard  alpha  distribution  [13]  as  the  initial  crack  size  distribution  were 
used  to  assess  the  MicroFaVa/DARWIN®  da/dN  and  life  predictions.  Essentially  identical 
results  were  obtained  by  the  current  methodology  and  those  reported  in  the  FAA  Circular.  In  the 
second  validation,  MicroFaVa  and  DARWIN®  were  used  to  analyze  the  probability  of  fatigue 
fracture  of  a  representative  Ti  rotor  design.  As  in  the  first  problem,  the  initial  crack  size  was 
described  by  the  distribution  of  the  hard  alpha  particle  size  and  a  power-law  FCG  relation  was 
used.  Figure  12(a)  shows  the  rotor  disk  design  divided  into  221  zones  after  5  zone  refinements. 
The  predicted  probability  of  fracture  computed  on  the  basis  of  conventional  (experimental  da/dN 
data)  material  input  and  MicroFaVa  input  are  compared  in  Figure  12(b).  It  should  be  noted  that 
the  results  based  on  the  conventional  data  input  did  not  consider  the  coefficient  of 
variation  (COV)  of  da/dN.  Without  considering  the  COV,  the  DARWIN®  prediction  using 
MicroFaVa  input  is  slightly  less  conservative  than  that  using  the  conventional  input.  The  ability 
of  MicroFaVa  to  predict  the  COV  in  addition  to  the  mean  da/dN  allows  one  to  quantify  the 
increase  in  the  probability  of  fracture  due  to  da/dN  variation.  As  shown  in  Figure  12(b),  the 
predicted  probability  of  fracture  is  increased  slightly  above  the  conventional  method  line  when 
the  COV  predicted  from  MicroFaVa  is  incorporated  into  the  DARWIN®  design  analysis.  These 
computations  illustrate  that  MicroFaVa  is  a  viable  methodology  for  predicting  fatigue  variability 
due  to  microstructural  variations  and  is  highly  compatible  with  probabilistic  design  and 
life-prediction  methods. 


14 


r — ii 

Refinement  5  -  221  Zones 

(a) 


Figure  12.  Analysis  of  a  Ti  rotor  design  using  MicroFaVa  and  DARWIN®:  (a)  hub 
geometry  with  the  critical  region  in  red,  and  (b)  computed  probability  of 
fracture  [16]. 
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4.0  CONCLUSIONS 


In  this  report,  the  development  of  a  probabilistic  microstructure-based  crack  initiation  and 
growth  model  was  summarized  and  the  model  capability  was  illustrated  for  life  prediction 
of  Ti-6A1-4V.  It  was  shown  that  the  microstructure-based  S-Nf  and  da/dN  simulation  results  are 
in  agreement  with  experimental  data  for  the  range  of  values  considered  in  this  study.  Fatigue  life 
variability  in  Ti-6A1-4V  appears  to  originate,  at  least  partly  if  not  totally,  from  microstructural 
variability  in  the  form  grain  size  variations.  In  addition,  the  microstructure-based  models  appear 
to  provide  reasonable  estimates  of  crack  growth  life  and  total  life  variability  when  compared 
with  life  values  based  on  experimental  da/dN  and  S-Nf  data.  Application  of  the  probabilistic 
fatigue  model  to  a  Ti  alloy  rotor  design  illustrates  that  the  proposed  model  is  capable  of 
predicting  the  confidence  bounds  of  the  fatigue  life  and  can  significantly  reduce  the  database 
requirement  used  to  define  fatigue  life  variability. 


5.0  LIST  OF  ACCOMPLISHMENTS 

The  accomplishments  achieved  in  this  program  are  as  follows: 

1.  The  development  of  a  set  of  physics-based  fatigue  models  for  predicting  fatigue  crack 
initiation  and  growth  in  structural  alloys.  This  suite  of  fatigue  models  allows 
predictions  of  FCG  rate  (da/dN  versus  AK)  response,  stress-life  (S-Nf)  data,  FCG 
thresholds  for  small  and  large  cracks,  and  dual  fatigue  limits  for  crack  initiation  and 
crack  growth. 

2.  The  formulation  of  a  probabilistic  framework  of  treating  fatigue  life  variability  due  to 
variations  of  microstructure.  This  formulation  enables  prediction  of  the  mean  and  the 
variance  of  fatigue  life  including  the  worst-case  limit  on  the  basis  of  the  probability 
density  distribution  of  the  pertinent  microstructural  length  scale  or  other  relevant 
material  variables. 

3.  Implementation  of  physics-based  fatigue  crack  initiation  and  growth  models  in  the 
probabilistic  framework  to  form  a  probabilistic  micromechanical  code  dubbed 
MicroFaVa  for  predicting  variability  in  initial  crack  size,  fatigue  crack  growth  rate, 
and  ultimately  fatigue  crack  initiation,  crack  growth  life,  as  well  as  the  total  life  and 
their  confidence  bounds  due  to  material  variability  or  uncertainties  in  material 
properties. 

4.  The  integration  of  MicroFaVa  with  a  probabilistic  design  and  life-prediction  code 
named  DARWIN®  (Damage  Assessment  of  Reliability  With  INspection)  and  the 
application  of  the  integrated  methodology  for  probabilistic  design  and  lifing  of  a  Ti 
rotor. 

5.  An  extensive  and  systematic  assessment  of  the  predictive  capability  of  the 
physics-based  fatigue  models  by  comparisons  against  a  large  experimental  database 
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of  fatigue  life  data  and  fatigue  crack  growth  rate  data  of  a  Ni-based  alloy  and 
a  Ti-based  alloy. 

6.  Verification  and  validation  of  the  integrated  probabilistic  design  and  lifing 
methodology  by  applying  the  computational  approach  to  two  benchmark  problems 
with  known  solutions.  The  utility  of  the  computed  material  input  has  been 
successfully  demonstrated  by  one  to  one  comparisons  against  the  conventional 
method  based  on  empirical  data  input. 

7.  The  application  of  the  probabilistic  micromechanical  code  to  elucidate  the  influence 
of  microstructural  variation  on  the  propagation  and  arrest  of  grain-sized  small  cracks, 
the  probability  density  function  of  small-crack  thresholds,  as  well  as  the  occurrence 
of  dual  fatigue  limits  for  crack  initiation  and  growth. 

8.  The  development  of  fundamental  understandings  on  several  microstructure/fatigue 
resistance  relationships:  (i)  the  dependence  of  the  intrinsic  threshold  on  the  planarity 
of  slip  and  dislocation  cell  formation,  (iii)  the  occurrence  of  dual  fatigue  thresholds 
and  dual  fatigue  limits,  (iii)  fatigue  life  variation  resulting  from  difference  in  the  slip 
morphology  and  enhanced  fatigue  crack  initiation  life  by  dislocation  cell  formation, 
and  (iv)  dependence  of  FCG  rate  on  dislocation  density  at  the  crack  tip. 
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3.  “Extension  of  a  Microstructure-Based  Fatigue  Crack  Growth  Model  for  Predicting 
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9.0  TRANSITIONS  AND  AWARDS 

9.1  Transitions 

The  technology  developed  in  this  program  has  been  incorporated  into  DARWIN  and 
transferred  to  the  USAF  DUST  program  (Contract  No.  F336 15-03-2-5203)  with  SwRI.  The 
Program  Manager  of  the  DUST  program  is  Dr.  Patrick  Golden  (phone:  (937)  255-5438)  at  Air 
Force  Research  Laboratory,  WPAFB,  Dayton,  OH. 

The  probabilistic  micromechanics-based  life  methodology  is  in  the  process  of  being 
transferred  to  a  DARPA-supported  program  entitled  “A  Revolutionary  Smart  Materials 
Approach  to  Structural  Health  Monitoring  and  Progress”  (Contract  No.  HR001 1-04-0004;  Dr. 
Leo  Christodoulu,  Program  Manager;  (703)  696-2374)  at  SwRI  and  an  AFOSR  program  entitled 
“Vehicle  Health  Monitoring”  (Contract  No.  FA955004 1 0254,  Dr.  Thomas  Kim,  Program 
Manager;  (703)  696-7312)  at  University  of  Texas  at  San  Antonio. 

9.2  Awards  Received 

No  award  was  received  during  the  reporting  period. 
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